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Abstract 

We compute the nonequilibrium real-time evolution of an 
0(A^)-~symmetric scalar quantum field theory from a systematic 
expansion of the 2PI effective action to next-to-leading order, 
which includes scattering and memory effects. In contrast to the 
standard expansion of the IPI effective action, the next-to- 

leading order expansion in presence of a possible expectation value 
for the composite operator leads to a bounded time evolution where 
the truncation error may be controlled by higher powers in 
We present a detailed comparison with the leading-order results 
and determine the range of validity of standard mean field type 
approximations . 

We investigate "quench" and "tsunami" initial conditions 
frequently used to mimic idealized far-from-equilibrium pion dynamics 
in the context of heavy- ion collisions. For spatially homogeneous 
initial conditions we find three generic regimes, characterized by an 
early-time exponential damping, a parametrically slow (power-law) 
behavior at intermediate times, and a late-time exponential approach 
to thermal equilibrium. The different time scales are obtained from a 
numerical solution of the time-reversal invariant equations in 1 + 1 
dimensions without further approximations. We discuss in detail 
the out-of-equilibrium behavior of the nontrivial n-point correlation 
functions as well as the evolution of a particle number distribution 
and inverse slope parameter. 

*Email: j.berges@thphys.uni-heidelberg.de 



1 Introduction and overview 

Current and upcoming heavy ion collision experiments at RHIC and the 
LHC have been an important motivation for the study of quantum field 
theory in and out of equilibrium. With these experiments one aims for 
a Quark Gluon Plasma produced in a transient state, finally releasing 
high multiplicities of particles of which the lightest hadrons are the pions. 
One major open question in the description of the underlying physics 
concerns the justification of current quantum field theoretical predictions 
based on equilibrium thermodynamics, local equilibrium or linear response 
assumptions. Immediately related to this appears the question for possible 
new signatures from processes which are far away from equilibrium. In 
contrast to a thermally equilibrated state, which keeps no information about 
the details of the initial conditions, nonequilibrium phenomena may help to 
understand in particular the earlier stages of a collision. 

A successful description of the dynamics of quantum fields away from 
equilibrium is tightly related to the basic problem of how macroscopic 
"dissipative" behavior arises from time-reversal invariant quantum dynamics. 
This is a fundamental question with most diverse applications. The 
techniques, which have to be developed to understand the physics of heavy- 
ion collision experiments, will be of relevance to the inflationary dynamics in 
the early universe, to the physics of baryogenesis or even to nonequilibrium 
aspects of mesoscopic quantum devices. 

The description of nonequilibrium quantum field theory from "first 
principles" can be based on a path integral formulation using the Schwinger- 
Keldysh techniques [|l|. In analogy to vacuum or thermal equilibrium 
quantum field theory one can construct a generating functional for 
nonequilibrium Green's functions which contains all quantum and statistical 
fiuctuations for given initial correlations or density matrix. 

However, controlled computational methods for the approximative 
solution of time evolution problems which respect the symmetries of the 
underlying theory are limited so far. In order to understand how macroscopic, 
effectively irreversible behavior arises one cannot use phenomenological 
approaches which typically break the time-reflection symmetry of the 
underlying theory. Simple (flnite-loop) perturbative descriptions can be 
unbounded and are known to break down at late times 0. While standard 
nonperturbative methods based on lattice Monte Carlo simulations have been 
successfully applied to thermal quantum fleld theory in Euclidean time, a 
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non-positive definite probability measure at real times already prevents their 
use. In contrast, classical field theory can be simulated. One only has to solve 
the classical field equations of motion sampled with the appropriate initial 
condition. Of course, the latter is not capable of treating genuine quantum 
effects but it should be reliable when the number of field quanta in each 
mode is sufficiently large. Another approach to which much work has been 
devoted restricts the discussion to systems close to thermal equilibrium or 
to effective descriptions based on a separation of scales in the weak coupling 
limit. For a recent review of the very interesting developments in this field 
seeRef. [|. 

Away from equilibrium, time-refiection invariant approximations to 
quantum field theories have been almost uniquely based on mean field 
type approximations (so-called "coUisionless" Hartree (-Fock) or leading 
order large- approximations) [Q. It is a known problem that the latter 
approximations are unable to describe thermalization at late time^, and 
their extensive use is based on the assumption that scattering effects may 
not change the picture dramatically for sufficiently early times. Efforts to go 
beyond mean field include a expansion of the generating functional for 
one-particle irreducible (1-P/) Green's functions beyond leading order 0. 
However, the standard expansion of the IPI effective action turns out 
to be unbounded in time [0]. For time evolution problems it therefore fails 
to give a controlled expansion. 

Recent progress came from a rather old concept based on a loop 
expansion of the generating functional for two-particle irreducible (2P/) 
Green's functions as introduced in Refs. [HI P, llOl O. It has been 



demonstrated in Ref. [|T^ for a scalar quantum field theory that a bounded 
time-evolution can be achieved without further approximations from a three- 
loop expansion of the 2PI effective action. The three-loop approximation 
reflects all symmetries of the underlying theory and goes beyond the mean 
fleld level, including scattering and memory effects. For very different 
nonequilibrium initial conditions it was shown from a numerical solution 
in 1 + 1 dimensions that the three-loop approximation leads to a universal 
asymptotic behavior of quantum flelds. The correlation functions at late 
time are uniquely determined by the initial (conserved) energy density 
and approach the equilibrium correlations of the corresponding three-loop 
thermal fleld theory [O]. 



^For improved results using inhomogeneous mean fields see Ref. 
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At three-loop order the 2PI effective action comprises the classical 
Boltzmann equation [|I4[ |T5|, |1^, 



16, 17, 1 



Without 



additional 

approximations the three-loop equations may be viewed as a "quantum 
Boltzmann equation" including off-shell effects and resumming an infinite 
order of derivatives [|l3l, which are difficult to include in kinetic descriptions 
T9| , pO| . However, despite the fact that the loop expansion reveals a 
consistent picture of the early- and late-time physics of quantum fields 
without further assumptions, in absence of a small expansion parameter it 
is not a controlled approximation. In this work we will discuss a controlled 
nonperturbative approximation. 



Our aim is to study a quantum field theory which captures aspects 
of low energy pion physics and which is simple enough that one can 
perform a quantitative treatment. The basic ingredient for the effective 
low energy description is the approximate chiral symmetry of the classical 
action for quantum chromodynamics (QCD). In QCD for Nj flavors the 
chiral SUl^Nj) x SUr^Nj) is broken to a vector-like SU{Nf), where the 
pions appear as the corresponding Goldstone bosons. For two flavors the 
SU{2) X SU{2) is locally isomorphic to an 0(4) symmetry which we will 
consider in the following. Unlike in chiral perturbation theory we will study 
here the corresponding linear model ("linear sigma model") for the light 
scalar and pseudo-scalar degrees of freedom which constitute the 0{N = 4)- 
vector {a, vf). For various this model has been extensively studied in 
thermal equilibrium, or nonequilibrium as in Refs. 0, ^0. 

We employ a systematic 1/N expansion of the 0(N)— symmetric 
2PI effective action r[0, G], which does not only depend on a possible 
expectation value of the quantum field, (p, but also on an expectation value 
of the time-ordered field bilinear, G (cf. Sects. ^ ^). Though = 4 may be 
used already as an expansion parameter, we will frequently consider larger 
values of A^ for the sake of a quantitative treatment. We compute the 
2PI effective action to next-to-leading order in 1/A^, which goes beyond 
the mean field-type leading-order approximation and includes scattering and 
memory effects. This approximation leads to a bounded time evolution for 
which the truncation error may be controlled by higher powers in 1/A^. In 
this case the employed expansion provides a controlled nonperturbative 
approximation which is not restricted to the weak coupling regime. Time- 

^For TV = 1 it also describes the long-range fluctuations of the QCD critical point pl[|. 
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Figure 1: For the 0(A^)-symmetric scalar theory in 1 + 1 dimensions we 
find three generic time regimes: An approximately exponential damping 
of oscillations characteristic for early times with inverse rate r^'^^'^^^ ~ 
X/^(damp)^ a parametrically slow or power law drifting at intermediate 
times, and a late-time exponential approach to thermal equilibrium with 
^(therm) ^ (therm) _ gQ^]^ for the employed "quench" and the "tsunami" 
initial condition t^'^'^^p) <^ ^(therm) None of the described different 
time scales are captured in the standard leading-order or mean field type 
approximations. The inclusion of scattering effects is crucial and the leading- 
order approximation is approximately valid only for t -C I/^^^'^'^p). 



reflection invariance and conservation of energy is preserved at each order 
in the expansion. For a discussion of similar approximation schemes 



applied to an A^-component quantum anharmonic oscillator see also Ref. |22 



For = 1 the current approximation includes the type of diagrams involved in 



the three- loop approximation of the 2PI effective action used in Refs. |12, 13 



The precise correspondence can be found by neglecting all contributions of 
the current approximation beyond three-loop and rescaling the four-point 



function "L L/3" in the three- loop approximation of Ref. |jT2 



We study the time evolution for two initial condition scenarios away from 
equilibrium. The flrst one corresponds to a "quench" (Sect. |7.1.1|) often 
employed to mimic the situation of a rapidly expanding hot initial state which 



cools on time scales much smaller than the relaxation time of the flelds |23|. 
Initially at high temperature we consider the relaxation processes following 
an instant "cooling" described by a sudden drop in the effective mass. The 
second scenario is characterized by initially densely populated modes in a 
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narrow momentum range around ±Pts (cf- Sect. [7.3. 1| ). This initial condition 
has been termed "tsunami" in Ref. p4| and is reminiscent of colhding wave 



packets moving with opposite and equal momentum. A similar nonthermal 
and radially symmetric distribution of highly populated modes may also be 
encountered in a "color glass condensate" at saturated gluon density with 



typical momentum scale pts ||2^) Of course, a sudden change in the two- 
point function of a previously equilibrated system or a peaked initial particle 
number distribution are general enough to exhibit characteristic properties 
of nonequilibrium dynamics for a large variety of physical situations. 

To reduce the numerical "cost" of following the evolution to sufficiently 
large times we present results for the 1 + 1 dimensional quantum field theory 
with spatially homogeneous correlations. We stress that in 1 + 1 dimensions 
many realistic questions can not be addressed, in particular since there is 
no spontaneous symmetry breaking. The 1 + 1 dimensional dynamics may 
help to understand the symmetric regime relevant for sufficiently high energy 
densities. However, the current results should be understood as a first 
controlled, quantitative approach to understand the dynamics of quantum 
fields far away from equilibrium. 

Independent of the details of the initial conditions we find the following 
characteristic time scales summarized in Fig. |l|. The nontrivial two- and four- 
point functions at next-to- leading order in 1/iV oscillate with initial frequency 
proportional to the renormalized initial mass Mjnit- All correlation functions 
quickly approach an exponentially damped behavior. A characteristic rate 
^(damp) obtained from the zero mode of the unequal-time two- 

point function F(t,0;p = 0) = (</)(t)0(O))p=o with a corresponding time 
scale r^*^^™P^ proportional to the inverse rate (Sect. |7.2| ). In this time 
range correlations with the initial state are effectively suppressed and 
asymptotically F{t, 0;p = 0) 0+. 

Both for the "quench" and the "tsunami" we observe that after the 
characteristic damping time the system is still far away from equilibrium. 
We find very different rates for damping and for the late-time exponential 
approach to thermal equilibrium]^, i.e. t^'^^'^p') <^ ^(therm) Though the 
exponential damping at early times is crucial for an effective loss of details 
of the initial conditions — a prerequisite for the approach to equilibrium 
— it does not determine the time scale for thermalization. At intermediate 



■^For the employed spatially homogeneous initial conditions we do not observe a power- 
law "tail" for the late-time evolution (see also Ref. pTj). 
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times the evolution of correlation functions is found to be characterized by 
a parametrically slow change (drifting) for times typically much smaller 
than We observe characteristic power law behavior \F{t,t;p) — 

F(0,0;p)| ~ i(p°™'''')p reminiscent of a hydrodynamic regime (Sect. |7.3|) . The 
large-time limit of correlation functions is determined by the (conserved) 
energy density and independent of the detailed initial conditions. 

Apart from the three qualitatively different time regimes, whose presence 
is rather insensitive to the details of the initial conditions, there are 
nonequilibrium phenomena specific to a "quench" or a "tsunami" (cf. 
Sect. 1^). After quenching from an initial thermal distribution with 
temperature Tq we find that the system relaxes asymptotically to an effective 
particle number distribution with temperature (inverse slope parameter) 
T < Tq. During the nonequilibrium evolution the nontrivial two- and four- 
point correlation modes show no substantially different behavior for low or 
for high momenta in absence of spontaneous symmetry breaking. In contrast, 
the "tsunami" initial condition can lead to strong differences in the dynamics 
for the low (p <^ pts) and the high lying modes (p ^ pts)- In particular, 
for the low momentum modes the inverse slope parameter quickly shoots 
up in response to the "tsunami" and can reach several Teq at intermediate 
times. The high momentum modes show the opposite effect and "cool" down 
slightly, before the inverse slope T(x°;p) approaches its asymptotic constant 
value Teq (cf. Sect. [7.4. 1|) . 

New phenomena of far-from-equilibrium physics include also the possible 
appearance of strongly interacting low momentum modes at intermediate 



times. For the employed "tsunami" of Sect. |7.3.1| the effective four-point 
function Xes{x^]P = 0)/QN reaches several times its initial value before it 
relaxes to its asymptotic late-time result, which is smaller that the initial 
value X/QN (cf. Sect. 7.3. 3| ). The strong renormalization of the effective 
four-point function at intermediate times also stresses the fact that effective 
descriptions based on a clear separation of scales can be difficult to achieve 
for far-from-equilibrium processes. 

Sects. and ^ present a short introduction to the generating functional 
for nonequilibrium 2PI Green's functions. In Sects. ^ and |^ we compute 
the 2PI effective action to next-to-leading order in and derive the 
equations of motion. Sect. |^ discusses the numerical implementation. The 
numerical results are given in Sect. |^, where we also present a comparison 
with leading-order results and determine the range of validity of mean field 
type approximations. We end with a conclusion and an outlook in Sect. ^. 
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2 2PI effective action 



We consider a quantum field theory with classical action S = J dfi^^x £, 

where ^ak^^-, x = (x°, x), is a real, scalar field with a = 1, . . . , N components 
(summation over repeated indices is implied). All correlation functions of 
the quantum theory can be obtained from the effective action r[0, G], i.e. 
the generating functional for two-particle irreducible (2PJ) Green's functions 
parametrized by the field Gab{x, y) representing an expectation value of the 
time ordered composite T{pa{x){pb{y) and the macroscopic field 0a (a;) given 
by the expectation value of ipa{x)- A discussion of the defining functional 
integral of the 2PI effective action can be found in Ref. |Q. Following ||^ it 
is convenient to parametrize the 2PI effective action as 

r[</),G'] = ^[(/)] + ^TrlnG'-^ + ^TrG'o'G' + r2[</),G'] + const (2.2) 

which expresses F in terms of the classical action S and correction terms 
including the function F2 which is discussed below. Here the classical inverse 
propagator iG^^^i^{x,y; (p) = 6'^S[(j)]/6(j)a{x)6(j)b{y) is given by 

^Go^abi^, y]4>) = - (^D,. + + A 0^(a;)0,(x)^ Sab^'^^^ {x - y) 

^ Ux)cPb{x)5''+\x - y) . (2.3) 



3iV 

In absence of external sources physical solutions require 

5F[0,G] 



50a(x) 

6Gab{x,y) 



0, (2.4) 
0. (2.5) 



Taking the derivative of ( p.2|) with respect to G one observes that the 
second stationarity condition ( p^.Sj ) is equivalent to the exact Schwinger- 
Dyson equation for the propagator: 

'^abi^^ y) = '^0 afe(^' ^) - ^ab{x, y, 0, G) (2.6) 
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with the proper self energy 

^ab{x,y;(f),G) = 2i—— -. (2.7) 

dGab{x,y) 

The simple relation between the self energy and the correction term r2 
motivates the representation ( p.2|) of the 2PI effective action. It is instructive 
to observe from ( \i.7\) that r2 must be two-particle irreducible (cf. also the 
detailed discussion in Ref. [§]). For ( |2.6| ) to be an identity S is to be 
interpreted as the proper self energy where only one-particle irreducible 
graphs can contribute. Suppose r2 has a two-particle reducible contribution 
of the form TGGT'. Then S has a contribution of the form TGT' since it 
is given by a derivative of r2 with respect to G. Such a structure cannot 
occur for the proper self energy and two-particle reducible contributions to 
r2 must be absent. 

Truncated Schwinger-Dyson equations are often used to describe 
nonperturbative physics. Typically the employed approximation is based on 
some ansatz for S which in turn corresponds to an efficient resummation of a 
large number of graphs. It is apparent from the simple relation ( [2.7] ) that the 
2PI effective action may be used as a powerful tool to obtain resummation 
schemes in a systematic way. 



The contribution r2 in ( |2.2| ) is given by all 2PI graphs with the 
propagator lines set equal to G P]. (A graph is two-particle irreducible 
if it does not become disconnected upon opening two propagator lines.) 
In presence of a nonvanishing expectation value the 2PI graphs 
are constructed from two kinds of vertices described by the effective 
interaction Lagrangian 

CwT{x;(f),<^) = -—(f)a{x)yDa{x)<^b{x)<^b{x) - ^(^(fa{x)<^a{x)^ ■ (2.8) 



The effective interaction can be obtained from (|2.1|) by shifting the field 
if ^ + (p. The terms cubic and higher in ip define the vertices from which 
the 2PI graphs contributing to T2 can be constructed. 

The relation of the 2PI effective action to the conventional generating 
functional for one-particle irreducible Green's functions is very simple. The 
IPI effective action corresponds to r[0, G] at that value of G for which ( |2.5|) 
holds, i.e. 

m-mc...] , mi„^,^.^,^o. (2.9) 



8 



3 Nonequilibrium Green's functions 



We will use the 2PI effective action to study the time evolution of n-point 
functions for given initial correlations at time tj. The initial n-point 
functions Tr{p(t/)(y9(t/, xi) . . .y9(t7,x„)} may include a mixed-state density 
matrix p{ti) for which TT{p^{ti)} < 1. This is necessary to allow for general 
nonequilibrium initial conditions and, in particular, if the initial density 
matrix describes thermal equilibrium. Often the initial conditions of an 
experiment — one may think of a scattering experiment of particles at 
an initial time tj far before they reach the interaction region — may be 
described by only a few lowest ra-point functions. For a setup which can be 
described in terms of an initial Gaussian density matrix the nonvanishing 
initial correlations are completely specified by and G and their first order 
time derivatives^ at tj, i.e. 

0(t/,x) = Tr{p{tiMtj,^)}, (3.1) 
0(t/,x) = Tr{p(t7)9^ov?(x°,x)}|^.o=f^ , 

and the initial two-point functions 

G(t,,x;t,,y) = TT{p{tiMtj,^Mtj,y)}-<f){ti,^)<f){ti,y), (3.2) 
H{ti,x;tj,y) = Ti{p{tj)d.^o(f{x°,x)Lp{tj,y)}l^o=tj - 0(t/, x)0(t/, y) , 
K{ti,x;ti,y) = TT{p{ti)d^oLp{x° ,x)dyO(^{y° ,y)}\^o=yO=tj - 0(t/, x)0(t/, y). 

We will assume in the following that the initial conditions can be described 
in terms of a Gaussian density matrix. It should be stressed at this point 
that a specific initial condition only restricts the experimental setup and 
represents no approximation for the time evolution. For t > tj irreducible 
higher n-point functions are induced by the time evolution of the interacting 
theory. The full 2PI effective action T[(j), G] contains all possible information 
about n-point functions for given initial conditions. 

A more general p(t/) involves the specification of nonzero higher initial 
n-point functions than (|3.1|) and (|3.2| ). The 2PI effective action in presence 
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For a Gaussian initial mixed-state density matrix there are three independent 



variances which can be chosen as G, H and K defined in (3^). Note that the 



initial correlation Tv{p{ti)(p{tj ,x.)dyOip{y^ ,y)}\y0^fj is not independent because of the 
commutation relation between ip and dyoLp. A detailed discussion of Gaussian initial 
conditions can be found in Ref. [ESl. 
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of a more general initial density matrix is discussed in |T0, O]. In general, 



a non-Gaussian p{tj) leads to additional source terms in the generating 
functional for 2PI Green's functions and F is no longer of the simple form 
( |2.2| ). These source terms are non- vanishing at time tj only and allow for the 
specification of the additional nonzero initial correlations. 

For Gaussian initial conditions the time evolution of n-point functions 
TY{p{ti)ip{x) . . . (p{y)} is completely specified by ( p.l|) and (|37^ ) and known 
F[0,G]. It should be noted that the n-point functions of the initial 
time problem differ from those used to compute S'-matrix elements in 
conventional scattering theory. The n-point functions for the initial time 
problem involve the trace over products of field operators. Stated differently, 
one computes matrix elements between same states. These are sometimes 
called "m — in" matrix elements in contrast to "m — out" matrix elements 
used in conventional scattering theory. The construction of the generating 
functionals for the respective n-point functions is identical up to the time 
path. The trace over products of field operators implies that the generating 
functional used for initial time problems involves a closed time path C 
The time path starts at the initial time tj and runs along the (real) time axis 
to some final time t and back to tj. The largest time of the path t is kept as 
a parameter and is evolved in the time evolution equations for the fields as 
described in Sect. 



4 expansion of the 2PI effective action 

In the interacting theory the 2PI effective action cannot be computed 
exactly. In this work we employ a expansion of the 2PI effective 
action F[0, G] where N is the number of field components. The leading order 
approximation is found by considering only those 2PI diagrams where the 
summation over internal lines yields a factor greater than or equal to A^' with 
/ the number of loops in the diagram. The next-to-leading order diagrams 
yield a factor N^~^, next-to- next-to-leading order A^'~^ and so forth. This 
defines a systematic classification in terms of powers of for 2PI diagrams 

^ Formulations of closed time path generating functionals often exhibit an infinite time 
interval ] — cx), oo[. In this case the number of field labels has to be doubled to distinguish 
the fields on the underlying closed contour. Causality implies that for any n-point function 
at finite times U {i = 1, . . . ,n) the contributions of the infinite time path employed in these 
formulations cancel for times larger than max(ti). For our purposes of a possible numerical 
implementation as an initial-value problem the infinite time path is not useful. 



10 




a b 



Figure 2: LO contribution to the 2PI effective action. 



contributing to r[0, G]. For a discussion of similar approximation scliemes 
applied to quantum mechanics see also Ref. p2 
be taken to be 0{^/N), hence Gab is 0{l) [ 



In the large- limit (pa can 
The leading order term is then 
proportional to and the next-to-leading order term is proportional to one, 
as will be considered in the following. 



4.1 2PI effective action to next-to-leading order 

For simplicity we will concentrate in the following on the symmetric regime 
where it is sufficient to consider r[0 = 0, G] = r[G]. We write 

r2[G] =r^«[G] + rfO[G] + ... (4.i) 

where denotes the leading order (LO) and r^'"^ the next-to-leading 
order (NLO) contributions. The dots indicate terms beyond the NLO 
approximation. The LO contribution to '^2[G\ is given by only one diagram 
shown in Fig. ||, 

Tf[G] = I d''+'xGaa{x,x)Gbb{x,x) . (4.2) 

Omitting combinatorial factors Fig. |^ indicates the diagrams contributing 
to r2 [G] in next-to-leading order. The first diagram of the series corresponds 
to the double bubble diagram shown in Fig. ^. It is followed by a three- 
loop diagram shown and then each subsequent diagram of the NLO series 
can be constructed from the preceding one by introducing at one vertex an 




Figure 3: NLO contribution to the 2PI effective action. 
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b a 

Figure 4: NLO "double bubble" contribution. 

additional loop as in Fig. ^. One may view these contributions as a three- 
loop diagram with an effective four-vertex containing a "chain" of bubbles. 
The double bubble and the infinite series of closed chain diagrams are the 
only graphs contributing to '^2[G\ to NLO. For instance, the other 2PI five- 
loop diagram shown in Fig. ^only contributes to NNLO. It is straightforward 
to sum the infinite number of NLO diagrams analytically. One finds 

r^LO[G'] = 1 /" d^+\ ln[B(G)](x,x) (4.3) 
2 Jc 

where 

B(x, y- G) = 5't\x -y)+^^ Gab{x, y)Ga,{x, y) . (4.4) 

It is instructive to expand the RHS of ([4.3|) to see the correspondence with 
the series of graphs in Fig. ^ 

j^d^+^x ln[B(G)](x,x) = j^d'^^^x (^t^Gab{x.x)Gab{x,x)) 

-- / d'^^^xd'^^^y{i—Gab{x,y)Gab{x,y)^{i—Ga'b\y,x)Ga'b\y,x)^ 

+ ... (4.5) 

The contributions ( |4.2| ) and ( |4.3| ) constitute the NLO approximation for 
T[G]. We will use this without further approximations to derive the time 
evolution equation for G in Sect. ^, which parametrizes the generating 



Figure 5: The five-loop 2PI "eye"-diagram contributes at NNLO. 
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functional T[G]. Once G is known for all times of interest the theory is 
"solved" in the symmetric regime to NLO: All correlation functions can be 
obtained from derivatives of the generating functional ( |2.2| ) and are given by 
(summed) diagrammatic expressions in terms of G. 



4.2 Comparison with conventional schemes 

We stress that in the above 2PI scheme the counting of orders of is 
done for diagrams with propagator lines associated to the field G. The time 
evolution for the propagator field is governed by a differential equation in G 
which is explicit in time and can be solved by standard techniques for given 



initial conditions (cf. Sect. |5.2|) . If not for a perturbative approximation there 
is no need to express G in terms of the classical propagator Gq. However, to 
compare with standard 1/N expansions for the IP I effective action G can 
be written as an infinite series parametrized by Gq. 

To leading order there is no difference to the conventional expansion 
scheme P]. The equivalence can be observed using ( |2.7D and iteratively 
expanding the result for the full propagator ( |2.6| ) in powers of the coupling 
A. This yields the solution of G in terms of an infinite series of diagrams with 
propagator lines associated to the classical propagator Gq. Replacing this 
solution for G in ( [4.2|) one recovers the conventional "bubble sum" expression 
for the leading order IPI effective action. 

Beyond leading order the counting for the 2PI effective action differs 
from the conventional expansion. The two schemes only differ by higher 
order terms. As an example, one may consider the three-loop diagram in 
Fig. H. From ( p.7|) one finds that this diagram leads to the "setting sun" type 
contribution to the self energy shown in Fig. |^, where lines correspond to the 
full propagator G. Expressing the resulting G in terms of Go one observes, 
in particular, the resummation of an infinite series of "ladder" diagrams 
where two lines repeatedly exchange momenta by one-loop subdiagrams 
connecting the two lines. This infinite series of diagrams is not included 
in the conventional expansion at any finite order in 1/iV. A detailed 
discussion of the importance of these diagrams in the context of transport 
coefficients can be found in Refs. 0, . 

In Sect. ^ we will observe that the employed approximation for the 2PI 
effective action leads to a bounded time evolution for the employed initial 
conditions, both at LO and NLO. In this case also the neglected term 
r2[G] - T^^[G] - F^LOfc] is finite at all times for a bounded full T[G] or 
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Figure 6: "Setting sun" type contribution to the self energy with hues G. 



T2[G\. The truncation error may therefore be controlled by its suppression 
with higher powers in This is in contrast to the conventional 1/N 

expansion for the IP/ effective action where terms grow unboundedly with 
time 0. In the case of contributions which grow proportional to powers of 
t"" /N for some a > such an expansion would break down at times of order 

The 2PI effective action is known to lead to a non-secular time evolution 
already for the three-loop approximation of the 2PI effective action as has 



been pointed out in Ref. [Q. In Ref. |22| the A^-component anharmonic 
oscillator, resumming an equivalent set of diagrams as discussed here, has 
also been shown to exhibit a bounded behavior. So far, we have not found 
initial conditions where this property could not be observed for the current 
approximation. To find the desirable analytical proof for arbitrary initial 
conditions is a difficult task in the presence of scattering and memory effects 
involved in the NLO time evolution, as is explained in the following. 

5 Equation of motion 

5.1 Schwinger— Dyson equation 

According to (|2.6| ) the exact inverse propagator can be written as 

I G-J{x, y) = -(□. + m^) 5ab 5't\x - - y- G) (5.1) 

where the full self energy S is given by ( p.7|) . In next-to-leading order in 
1/A^ the approximation for the self energy follows from the derivative of 
(113), O with respect to Gab- One finds 

E„i,(x, y] G) = G^c{x, x)5ab Sc'^^{x - y) - Gab{x, y)B-'^{x, y; G) 

(5.2) 

with B given by (^4.4|) . Multiplying (^4.4|) with B ^ and using the identity 

j^d'''^^zB{x, z)B^^{z,y) = 5'^~^^{x — y) one observes that the inverse of B 
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obeys 

B-\x,y;G) = 5t\x - y) - i I{x,y) . (5.3) 

Here the function /(x, y) resums the "chain" of bubble graphs discussed in 
Sect. i.e. 

A A /" 

I{x^y) = -^Gabix,y)Gabix,y) - i— d'^^'^z I{x,z)Gab{z,y)Gabiz,y) . 
bJM bJM Jc 

(5.4) 

Rewriting ( p.3[ ) and ( |5.4| ) we note that the combination 

^ B-\x, y) = -^ {^'-''i^-y)-^-^ J d'^'z B-i(x, z)Gab{z, y)Gab{z, y)) 

(5.5) 

plays the role of the (nonlocal) four-point function to order 

In the symmetric regime we can evaluate ( ^.1|) for the configuration 
Gab{x,y) = G{x,y)6ab- In this form the C(l) and 0{1/N) contributions 
are explicit. We separate E in a local part and a nonlocal part, 

E{x,y;G) = -iS('°="i)(a;; G')5^+^(x - y) + S(°°°'°^"i)(a;, G) . (5.6) 

Since corresponds to a local mass shift it is convenient to write 

M\x; G)=m^ + j:^^°^^^\x; G). (5.7) 

With this notation the Schwinger-Dyson equation for the propagator reads 

iG-\x, y) = - (a + M\x- G)) 5^+1 (x - y) - i S(-'^i°^-')(x, y- G) . (5.8) 

From ( [5.2| )-( p7^ ) we observe that the local part receives LO and NLO 
contributions, 

M\x- G)=m' + A^^ G(x, x) , (5.9) 
while the inhomogeneous part of the self energy is nonvanishing only at NLO, 

S(-°'°-i)(a:,y;G) = - A g(x, i/)/(a:, y) , (5.10) 

where the resummed "chain" of bubble graphs ( |5.4| ) reads 

I{x,y) = ^G{x,y)G{x,y)-i^ [ d^+^z I{x,z)G{z,y)G{z,y) . (5.11) 



It is important to note that for the LO, or ^ oo, approximation only 
a local mass shift is induced by the interactions. The NLO contribution is 
crucial in order to include direct scattering. 
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5.2 Initial- value time evolution 

The form of the equation of motion ( |5.8|) is suitable for a boundary value 
problem and can be used, in particular, to discuss the propagator in 
thermal equilibrium. However, nonequilibrium time evolution is an initial 
value problem. Here a suitable form of the equation of motion can be 
obtained by multiplying ( |5.8D by G from the right and by integration. With 
df^+^z G-\x, z)G{z, y) = 6^-^\x - y) one finds 

+ M2(x; G))G{x,y)+i J d'^^hj:^''°''^°'''^\x, z;G)G{z,y) = -i5'^+\x-y) 

(5.12) 

This represents a second order differential equation for G which may be used 
to study the time evolution for given initial conditions ( |3.2| ). Very useful 
further simplification arises if we decompose 

G{x,y) = G>(x,y)ec(a;°-/)+G<(x,y)ec(2/°-x°), 
j.(noniocai)(^^^) = S> (x, y) Oc (a;° - + S< (x, y) Be - 0;°) (5.13) 
and observe 

n^G{x, y) = ec(xO - /)a.G>(a;, y) + Qciv'' - a;°)aG'<(x, y) - z5^+i(x - y) . 

(5.14) 

Note that the equal-time commutation relation for the fields indeed imply 

^ (g> (x, y) - G< (x, y)) ^^^^^^^ = -15'' (x - y ) . (5.15) 

For a real, scalar field theory Gy{x,y) = G<{x,y) = Gy{y,x) and all 
equations can be expressed in terms of one function (G>) only. The self 
energy has the same properties S>(x, y) = S<(x,y) = S>(?/,x), which can 
be verified from ( [5.10| ) and is valid to all orders in 1/A^. The 6-functions in 



(|5.14| ) are defined along the time contour G and can be evaluated 

j dz^T.y{x,z-G)G>{y,z) 



+ j dz^J:y{x,z;G)Gy{z,y) - j dz^ J:y{z,x;G)Gy{z,y) 
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This yields the evolution equation 

{n^ + M\x-G))G^{x,y) = 2 j t/zj (5.16) 

x° y° 

j dzHm[J:>{x,z;G)]G>{z,y) - j dz^ J:>{x, z;G) lm[G>{z,y)] 



We note that ( ^.16] ) represents an exact equation for known local part 



M^(x; G) and nonlocal part S>(x, ?/; G) of the self energy 
5.3 Real-valued evolution equations 

For the real scalar field theory there are two independent real-valued two- 
point functions, which can be associated to the expectation values of the 
commutator and the anti-commutator of two fields. We note that the field 
anti-commutator is determined by the real part and the commutator by the 
imaginary part of the complex function G^{x,y). Following Ref. we 
define 

F{x,y) = ^(GUx,y) + Gl{x,y)^=Re[G^{x,y)], (5.17) 

p{x,y) = i(Gy{x,y)-Gl{x,y)^=-2lm[G^{x,y)] (5.18) 

where F is the "symmetric" propagator and p denotes the spectral function, 
with the property F*{x,y) = F{x,y) = F{y,x) and p*{x,y) = p{x,y) = 
—p{y,x). Similarly, we define 

^F{x,y) = ^(s>(x,i/) + S;(x,i/)) =Re[S>(x,y)], (5.19) 

Sp(x,y) = «(s>(x,i/)-S;(x,y)) =-2Im[S>(x,i/)]. (5.20) 
The time evolution equations for F and p follow from ( |5.16|) |T3[ 

+ M\x- F)) F{x, y) = -Jdz!^ (5.21) 

xO y° 

j dz^ T.p{x,z;F,p) F{z,y) - j dz^ T.f{x, z; F, p) p{z,y)^ , 
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(a + M^x; F)) p(x, y) = -jdzj dz' z; F, p) p{z, y) . (5.22) 

yO 

The form of the evolution equations for the symmetric propagator and the 
spectral function is independent of the approximation. They are equivalent 
to the exact Schwinger-Dyson equation ( |5.1| ), however, in the form ( |5.21| ) 
and (|5.22|) they are more suitable for initial time problems. 

We note that due to the canonical commutation relations or ( |5.15|) the 
spectral function obeys the equal-time properties 

p{x,y)\^o=yO = 0, d^,op{x,y)\^,o=yO = 5'^(x-y). (5.23) 

The time evolution of the spectral function is completely determined by (|5.22|) 
and (lOl . 

In the NLO approximation M^(x; G) is given by (|5.9|) , i.e. 



M\x;F) = m^ + X^^^4r^F{x,x) , (5.24) 
6N 



and from ( |5.10|) and (|5.11|) one finds 

A /„ 1 



^F{x,y;F,p) = -— (^F{x,y)lF{x,y) - -p{x,y)Ip{x,y)Y (5.25) 
^p{x,y;F,p) = —^(^F{x,y)Ip{x,y)+p{x,y)lF{x,y)^ (5.26) 

with^ 



lFix,y) = ^(^F\x,y)-^p'ix,y)) (5.27) 

I dz'lpix,z)[F\z,y)-^p'{z,yi 







-2 J dz'lF{x,z)F{z,y)p{z,y) 







^Note that the resummed "cham" of bubbles ( 5.11 ) can be written as I{x,y) = 
lUx,y)ecix° - 2/0) + I<ix,y)eciy° - x") with r>(a:,y) = J<(a:,y) = I>iy,x). With 
this notation lF{x,y) = ^(I>{x,y) + /*(a;,y)) = Re[/>(a;,y)] and Ip{x,y) = i{Iy{x,y) - 
/*(x,y)) = -2Im[/>(a;,y)]. 
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Ip{x,y) = ^F{x,y)p{x,y) (5.28) 



X" 



We emphasize that F, p,TiF and Sp are real functions, and one observes that 
the evolution equations ( p.21[ )-( p.28[ ) are causal and time-reversal invariant. 



5.4 Comparison with thermal equilibrium 

It is instructive to consider for a moment the quantities F and p, S^? 
and Ep, as well as Ip and Jp in thermal equilibrium. The 2PI effective 
action in thermal equilibrium, T^'^'^\ is given by the same expression ( p.2|) , 
with ( [4.2|) and ( |4.3|) in the NLO large-A^ approximation, if the closed time 
path is replaced by an imaginary path C = [0,—i(3]. Here (3 denotes the 
inverse temperature. Since in equilibrium these functions depend only on 
the relative coordinates it is convenient to consider the Fourier transforms 
F(^<^)(^,p), p(^^)(^,p), 4'^(a;,p), S^)(.;,p), I^P\u;,p) and 4^^)(^,p). 
From the periodicity ("KMS") condition for the propagator in imaginary 
time, G{x,y)\xO=o = G{x,y)\x0=-ii3, one infers the generic equilibrium 
relation^ (cf. also Ref. 0) 



F^^^\u,p) = -^(n^{uJ) + ^-)p'^^^\uJ,p), (5.29) 
E^P\co,p) = -t{n^{co) + ^)j:f^\co,p), (5.30) 
I^P\u;,p) = -^(n^(u;) + ^)Ii:^^\u,p) (5.31) 



with nB{uj) = (e^'^ — 1) ^. Here we have used that with (^.11|) and (|5.1(]| ) the 



functions /(x, y) and xi('^°'^^°'^^')(x, y) satisfy the same periodicity condition as 
G{x,y). 

While the spectral function p'-^'^^ encodes the information about the 
spectrum of the model, one observes that the symmetric function F^^"^^ 
encodes the statistical aspects in terms of the particle distribution function 
ub and similarly for the E^'^'^^- and J ^^'^^ -functions. We note that the ratio 

^In our conventions the Fourier transforms of the real-valued antisymmetric functions 
p{x,y), Tip{x,y) and Ip{x,y) are purely imaginary. 
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Ep'"'^(a;, p)/2a; plays in the limit of a vanishing w-dependence the role of a 
decay rate for one-particle excited states with momentum p. The function 
Ip'^'^\uj, p) encodes information about the effective vertex as will be discussed 
in Sect. [7.2.3| . In the following we return to the nonequilibrium case and study 



the time evolution equations of Sect. 



6 Numerical implementation 



The time evolution equations (|5.21| )- (|5.28| ) are nonlinear, integro-differential 



equations. Though the equations are in general too complicated to be tackled 
analytically, they can be efficiently implemented and solved on a computer. 
Here it is important to note that all equations are explicit in time, i.e. all 
quantities at some later time tf can be obtained by integration over the 
explicitly known functions for times t < tf with given initial conditions. This 
aspect is discussed in more detail below. In this respect, solving the initial 
time evolution equations for nonequilibrium problems is simpler than solving 
the corresponding gap equation of the form ( |5.8|) self-consistently, which is 
typically employed for the study of thermal equilibrium. For simplicity we 
consider spatially homogeneous fields where 

y) = F{x', j/O; X - y) = 1 ^ e^^^^'y^ F{x' , y'; p) (6.1) 

and similarly for p{x,y) = p{x^,y^;:si — y). Spatially inhomogeneous fields 
pose no additional complication but are computationally more expensive. 

We employ a time discretization x° = nat, y^ = mat with stepsize at, 
F(x°, y^; p) = F{n, m; p) and 

dloF{x^,y^) ^ \(F{n + l,m) + F{n-l,m) -2F{n,m)), (6.2) 
at \ / 

/n—1 
dtF{t,y°) ^ at(^F(0,m)/2 + ^F(Z,m) +F(n,m)/2). (6.3) 

'=1 

The above simple discretization leads already to stable numerics for small 
enough stepsize at (cf. also Refs. |0 |1^), but the convergence properties 
may be easily improved with more sophisticated standard estimators. 

The discretized equations for the time evolution of the matrices F{n,m) 
and p{n, m) advance time stepwise in the "n-direction" for fixed m. As for the 
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continuum the propagators obey the symmetry properties F{n, m) = F{m, n) 
and p{n,m) = —p{m^n). Consequently, only "half" of the (n, m)-matrices 
have to be computed and p(n, n) = 0. Similarly, since Ip is antisymmetric in 
time Ip{n, n) vanishes identically. It is also useful to note that p(n+l, n) = at 
is fixed by commutation relations (cf. Eq. (|5.23|) ). As initial conditions 
one has to specify F(0,0;p), F(l,0;p) and F(l,l;p), while p(0,0;p) = 
p(l, 1; p) =0 and p(l, 0; p) = at are fixed. 

The time discretized versions of the evolution equations (|5.21|) and ( ^.22|) , 

F{n + 1, m; p) = 2F{n, m; p) — F{n — 1, m; p) 



2)2, 2 , ^ ^ + 2 



-a; ip-" + + X 



F{n, n; k) j F(n, m; p) 



-at 



6N 

S,(n, 0; p) F(0, m; p)/2 - S^(n, 0; p) p(0, m; p)/2 (6.4) 



m— 1 



+ Yl (^p(^' ^' P) ^(^' P) ~ ^F{n, I; p) p(/, m; p; 
1=1 

n— 1 s 

+Y1 ^p(^'^'P)^(^'"^'P) f ' 

l=m ^ 

(ra > m) and similarly for p are explicit: Starting with n = 1, for the time 
step n + 1 one computes successively all entries with m = 0, . . . , n + 1 from 
known functions at earlier times. At first sight this property is less obvious 
for the non-derivative expressions ( |5.28| ) for Ip and ( p.28| ) for Jp whose form 
is reminiscent of a gap equation. However, the discretized equation for Jp, 

N ^ M'^k r , , , 

/p(n, m; q) = - y |^(^' Q - k)p(n, m; k) 

n— 1 X 

^ /p(ra,/;q)F(/,m;q-k)p(/,m;k) > , (6.5) 

— ^_i_i 



-a* 

/=m+l 



shows that all expressions for Ip{n,m) are explicit: Starting with m = n 
where Ip vanishes one should lower m = n, . . . ,0 successively. For m = n — 1 

*For the discretization of the time integrals it is useful to distinguish the cases n > m 
and n < m. We compute the entries F{n+ 1, m) from the discretized equations for n > m 
except for n + 1 = TO where we have to use the equations for n < m. 
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one obtains an explicit expression in terms of F{n, m) and p(n, m) known 
from the previous time step in n. For m = n — 2 the RHS then depends on 
the known function Ip{n,n — 1) and so on. Similarly, for given Ip[n,m) it 
is easy to convince oneself that the discretized Eq. ( |5.28| ) specifies lF{n,m) 
completely in terms of 0), . . . , Ip^n, m—1), which constitutes an explicit 
set of equations by increasing m successively. 

It is crucial for an efficient numerical implementation that each step 
forward in time does not involve the solution of a self consistent or gap 
equation. This is manifest in the above discretization. The main numerical 
limitation of the approach is set by the time integrals ( "memory integrals" ) 
which grow with time and therefore slow down the numerical evaluation. 
Typically, the influence of early times on the late time behavior is suppressed 
and may be neglected numerically in a controlled way pO|. 



6.1 Continuum and infinite volume limit 

We use a standard lattice discretization for a spatial volume V with periodic 
boundary conditions. For a spatial volume V = [NgasY with lattice spacing 
ttg and periodic boundary conditions momenta are given by 



d 



2 4 . 2 fasPi\ l-nui 



=1 



s 



al V 2 / ' ^' iV,a 



where = 0, . . . , A^s — 1. The lattice introduces a momentum cutoff vr/a^. 
In order to obtain the theory in the continuum the lattice spacing has 
to be sent to zero and renormalization is necessary. In the following we 
consider numerical results in 1 + 1 dimensions. In this case the employed 
renormalization condition for the mass term, M^(a;) l^^o^o = j^init; 
described by ( |5.24| ) is sufficient to obtain finite, physical results. 



In order to study the infinite volume limit one has to remove finite 
size effects. Here this is done by increasing the volume until convergence 
of the results is observed. Fig. shows a typical time evolution of the 
unequal-time correlation zero mode 0;p = 0) for three different volumes 
L* = iVsOs = 256 X 0.4/MiNiT, -Z^*/4 and L*/16. The initial conditions are 
the same as for Fig. |10| described in detail below. One observes that the 
damping tends to be stronger for smaller volumes. We find that increasing 
the volume beyond L^, does not change the results such that one could 
distinguish the curves in the plot (< 1%). Depending on the initial condition 
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Figure 7: Volume dependence of the effective damping of oscillations, shown 
for the zero mode of the unequal-time function F{x^,0;p) in the NLO 
approximation for three different volumes. Here = 256 x 0.4/Minit 
denotes the "infinite" volume limit. The initial conditions are the same as 
for Fig. and the evolution is discussed in detail below. 



and the considered quantity the employed volumes can, however, be as large 
as L = 1600 X 0.4/MiNiT as for the LO evolution shown in Fig. |^, or as small 
as L = 32 X 0.2/MiNiT employed in the second frame of Fig. |13| showing the 
large-time behavior of the effective four-point function. (For the first frame of 
Fig. showing the coupling at early times, we use L = 256 x 0.2/Minit and 
one observes a small volume dependence for the large-time behavior of this 
quantity.) We vary our lattice spacings in the range a<j = (0.2 — 0.4)/Minit 
and at /as = 0.05 — 0.2.0 

^ For time evolution problems the volume which is necessary to reach the infinite 
volume limit to a given accuracy can depend on the time scale. This is, in particular, 
due to the fact that finite systems can show characteristic recurrence times after which 
an initial effective damping of oscillations can be reversed. The observed damping can 
be viewed as the result of a superposition of oscillatory functions with differing phases or 
frequencies. The recurrence time is given by the time after which the phase information 
contained in the initial oscillations is recovered. Then the damping starts again until 
twice the recurrence time is reached and so on. In 1+1 dimensions we explicitly verify 
that in the LO approximation {N oo) the observed recurrence times for the equal-time 
correlation F(x, x) scales with the volume L = N^Gs or the number of degrees of freedom 
to infinity. We emphasize that the phenomenon of complete recurrences, repeating the full 
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7 Nonequilibrium dynamics: damping, drifting and 
thermalization 



The spatially homogeneous solution of ( |5.21| ) or ( p.22| ) in the limit of a 



free field theory (A = 0) describes modes which oscillate with frequency 
a/p^ + w? /2t: for unequal-time functions F(a;°,0;p) and p(a;°,0;p), unless 
they are not identically zero. Equal-time correlation modes F(x*^,x°;p) 
oscillate either with twice that frequency or they are constant in timeQ. The 
latter corresponds to solutions which are translationally invariant in time, 
i.e. ?/°; p) = — ?/°; p) and ?/°; p) = p{x^ — y^; p). 

Time translation invariant solutions play an important role for the 
dynamics of nonequilibrium field theory. In the presence of interactions we 
will observe below that for a large variety of initial conditions all modes 
approach time translation invariant solutions at asymptotically large times. 
This qualitative property of the asymptotic solution can be observed both for 
the LO and the NLO approximation. However, at LO the large-time behavior 
depends explicitly on the (conserved) initial particle number distribution and 
thermalization cannot be observed as is shown below. In contrast, in the 
NLO approximation scattering is taken into account and the nonequilibrium 
late time result turns out to be insensitive to the detailed initial conditions, 
approaching a thermal equilibrium distribution which is determined by the 



(conserved) initial energy density. Following Ref. we will call in the 
following time translation invariant solutions of the evolution equations fixed 
point solutions. 

We stress that during the nonequilibrium time evolution there is no loss 
of details about the initial conditions in any strict sense. The evolution 
equations are time reflection invariant and at any time the evolution can 
be reversed and the initial conditions recovered]^. In particular, starting 
away from a time translation invariant solution no fixed point solution can 
ever be reached during the nonequilibrium evolution (cf. also the discussion 



initial oscillation pattern after some characteristic time, is not observed once scattering 
is taken into account. Periodic recurrences can occur with smaller amplitudes as time 
proceeds and are effectively suppressed in the large-time limit even for small systems. As 
a consequence, in the NLO approximation we find that convergence of results to a given 
accuracy can be obtained for sufficiently large, fixed volumes ind epen dently of the time. 



Note that for equal times p{x^,x^;p) = according to Eq. ( 5.23 ). 
In this paper we consider closed systems without coupling to an external heat bath 



or external fields, which could provide sources or sinks of energy. 
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in Sect. |7.1.TD . However, fixed point solutions can be approached arbitrarily 
closely for sufficiently large times as is discussed below. 

In the following we present our numerical results and provide a detailed 
explanation of the summary in Sect. |l]. To make contact with earlier work 0] 
we begin with a short discussion of the dynamics for the LO approximation 



which can be characterized by a fixed point solution |32]. The LO fixed 
point is unstable once scattering is taken into account and we present the 
comparison with the NLO approximation afterwards. 



7.1 The LO fixed point 

The LO contribution to the effective action ( [4.2|) adds a time dependent 
mass shift to the free ffeld evolution equation. The resulting effective mass 
term M^(x), given by ( ^.24|) for oo, is the same for all Fourier modes 

and each mode propagates "coUisionlessly" . The evolution equation (|5.21[) 



for F{x,y) becomes nonlinear. However, it remains local in time since the 
"scattering" contributions on the RHS are absent in the LO approximation. 
In this case one observes that the evolution of F decouples from p. In LO 
the spectral function does not inffuence the time evolution of the symmetric 
propagator similar to the free field theory limit. 



7. LI "Quench" 

In Fig. H we plot the evolution of M^(x) in the LO approximation as a 
function of time The time evolution follows a "quench" described by an 
instant drop in the effective mass term from = 2M^^^r^ to M^j^^r^. (Cf. 
the discussion in Sect. ^) The initial equilibrium high temperature (Tq = 
2Minit) particle number distribution is no{p) = l/{exp[^/Jp^f^^^ /Tq] — 1). 
The sudden change in the effective mass term drives the system out of 
equilibrium and one can study its relaxation. We present M^(x) for three 
different couplings A = Aq = O.SMpj^jrp (bottom), IOAq (middle) and 40Ao 
(top). For the plot all quantities are given in units of appropriate powers 
of the initial time mass Minit = M{x)\xO=q. One observes that M'^[x) 
shoots up from Mi^^r^ in response to the "quench". The initial oscillations 
with frequency Minit/tt are damped and at large times M^(x) approaches a 
constant M|p depending on tiq and A. 

If a strictly time translation invariant solution were reached, with 
F{x,x) = F{x — x), then the LO mass term would be a constant and 
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Figure 8: Shown is the mass term M'^{x) in the LO approximation for 
three different couphngs A = Ao,10Ao,40Ao in units of appropriate powers 
of Mj^j^jrp. The asymptotic constant M|p depends on the couphng and the 
initial particle number distribution uq. The value Mpp is well described by 



the LO gap equation (|7. 1|) . The nonthermal LO fixed point becomes unstable 
once scattering is taken into account in the NLO approximation. 



given by = + ^F{x - x = 0) = Mj^it according to Eq. (|5.24|) 



for — ^ oo. Clearly for the initial conditions employed in Fig. ^ the 
mass term is non-constant and the asymptotic values deviate from Mi^irj,. 
We may compare the asymptotic dynamics with the solution of the LO 
equations for given particle number distribution no(p) and an initial mass 
term such that M'^{x) = Mq^p is constant. For a constant mass term 
the LO approximation behaves similar to a free theory. The propagator 
modes — y^]p) = [no{p) + 1/2] /lj{p) cos [uj{p) (x° — ?/°)] oscillate with 
uj = a/p^ + Mq^p, where the effective mass term according to Eq. ( [5.24} ) is 
given by the LO gap equation^ 



r2 2 , ^ f dp ( , , , 1 \ 1 



Here the logarithmic divergence of the one- dimensional integral is absorbed into the 
bare mass parameter m? using the same renormalization condition as for the dynamical 
evolution in the LO approximation, i.e. ™^ + f / ('^o(p) + \) [p^ + Mq )^^/^ = M'^^yt ■ 
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The result from this gap equation is M^p^-p = {1.01, 1.10, 1.29}Mpj^jrp for the 
three values of A. For this wide range of couplings the values are in good 
numerical agreement with the corresponding dynamical large time results 
inferred from Fig. | as M|p = {1.01, 1.11, 1.31}MiV at x° = 200/Minit. 
We explicitly checked that at = 400 these values are still the same. 

We stress that the particle number does not change during the LO time 
evolution and no approach to equilibrium is observed. Indeed the particle 
number 

nip) + l = (F(x°,i/°;p)a,o9,oF(xO,y°;p)-(9,oF(x°,i/°;p))')\o=,o (7.2) 

is a strictly conserved quantity in the LO approximation ^, For the 
quench discussed above dxoF{x^,y^]p)\xQ=yO = with n{p) = no{p) and the 
initial high temperature distribution remains constant. 

We conclude that for the spatially homogeneous situation considered here 
the LO approximation cannot describe the approach to thermal equilibrium. 
The late time behavior is characterized by the LO fixed point which explicitly 
depends on the initial particle number distribution. This nonthermal fixed 
point becomes unstable once scattering is taken into account as is discussed 
in the following. Before we address the nonequilibrium dynamics in the NLO 
approximation in more detail we note from Eqs. ( [5.21| ) and ( |5.22| ) that the 
LO approximation for F{x,y) and p{x,y) becomes exact^ for x^,y^ —>■ 0. 
For very early times one therefore expects the LO approximation to yield a 
quantitatively valid description. This expectation is the basis for the many 
applications of the LO approximation found in the literature so far. Within 
the present framework we can test when the LO approximation breaks down 
by comparing it with the NLO results. We will see that the LO approximation 
fails to describe the main qualitative aspects of an early time exponential 
damping or late time approach to thermal equilibrium. However, for weak 
enough couplings certain time-averaged quantities at early times can be well 
determined by the LO fixed point solution. The particular role of the LO 
fixed point solution for time-averaged quantities at early times was pointed 



out for the classical field theory limit in Ref. 32 



This is due to the fact that we choose to consider Gaussian initial conditions, where 
at a;*^ , y*^ = irreducible n-point functions vanish identically for n > 2. 
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Figure 9: Time dependence of the equal-time zero mode = 0) 

after a "quench". The employed coupling is X/6N = {5/6N ~ 0.083) M^^-^rj, 
for = 10. The inset shows F{x, x) which is the sum over all momentum 
modes. (All in units of appropriate powers of Minit.) 



7.2 Early-time exponential damping 

In contrast to the LO approximation at NLO in the expansion the 
evolution equations ( |5.21 ) and (|5.22 ) include the nonlocal parts of the self 



energy, 'Ep and S^, which encode scattering effects. In particular, the 
equations become nonlocal in time and the evolution depends on the time 
history. As pointed out above these memory effects can already be observed 
on the level of the exact equations ( |5.21|) and (|5.22|) for known T,F,p- The time 



structure of the equations reflects causality and time reflection invariance 
which is preserved in the expansion scheme. 

In Fig. ^ the solid line shows the behavior of the zero momentum mode 
of the equal-time correlator F{x^,x^;p) in the NLO approximation after a 
"quench". The inset displays the sum over all equal-time modes, F{x,x) = 
J and M^{x)\^o=o = ^iV- The coupling is X/QN = {5/6N ~ 

0.083) Mpj^jrp for = 10 and the initial particle number distribution no(g) is 
the same as described in Sect. [7.1| for the LO approximation. For comparison 
the dotted lines in Fig. show the corresponding LO evolution. 
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Figure 10: Shown is the evolution of the unequal-time correlation 
F(a;'', 0; p = 0) after a "quench". Unequal-time correlation functions 
approach zero in the NLO approximation and correlations with early times 
are effectively suppressed. (In units of Minit-) 

7.2.1 Damping 

As expected from the previous discussion, for very early times the NLO 
behavior is well described by the LO approximation. However, one observes 
that the damping of oscillations is more efficient at NLO due to scattering 
effects. The amplitude of the zero mode F{x^,x^;p = 0) in Fig. ^ at NLO 
shows a substantial difference to the LO one. Summed over all momentum 
modes, however, the difference becomes less as can be inferred from F{x, x) 
shown in the inset of the Figure. The different damping behavior becomes 
even more pronounced for unequal-time correlators. In Fig. |1^ we plot 
F{x^,0]p = 0) as a function of time. In sharp contrast to the NLO result, 
where the maximum amplitude of the correlation mode is damped, one 
observes that the LO correlator is not damped at all — in accordance with the 
LO fixed point behavior described in the previous section. We find that the 
unequal-time two-point functions approach zero in the NLO approximation 
and correlations with early times are therefore suppressed once scattering 
is taken into account. We stress that time- reversal invariance implies that 
the oscillations can never be damped out to zero completely during the 
nonequilibrium time evolution, however, zero can be approached arbitrarily 
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Figure 11: The correlation modes oscillations quickly approach an 
exponentially damped behavior. The logarithmic plot of \p{x^,0;p = 0)| 
and \F{x^,0]p = 0)| as a function of time x° shows an oscillation envelope 
approaching a straight line (in units of Minit). 



closely. 

In the NLO approximation we find that all modes F{x^,y'^]p) and 
p{x^,y^]p) approach an approximately exponential damping behavior for 
both equal-time and unequal-time correlations. In Fig. |TT| the approach 
to an exponential behavior is demonstrated for the unequal-time functions 
\p{x^,0]p = 0)1 and \F{x^,0]p = 0)|. The logarithmic plot shows that after 
a non-exponential period at very early times the envelope of oscillations can 
be well approximated by a straight line. From an asymptotic envelope fit of 
F{x^, 0;p = 0) to an exponential form ~ exp(— 70'^'^ t) we obtain a damping 

rate 7^'^'"^^^ = O.OIGMinit- 

We conclude that, though strict dissipation can never be observed in a 
time-reversal invariant and energy conserving system, characteristic features 
of dissipative systems like exponential damping can be realized to very 
good approximation. The fact that correlation functions show an effective 
exponential damping once scattering is taken into account seems not to 
depend on the details of the approximation. We have observed the same 
qualitative feature in a three-loop expansion of the 2PI effective action 
0,0. 
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7.2.2 Effective mass 

One observes from Fig. ^ that the correlation mode F{x^,x^]p = 0) shoots 
up after the "quench" and relaxes. The initial oscillation frequency of the 
equal-time mode, Minit/tt, changes only slowly. The early-time behavior 
of the NLO mass term M^(x) is rather well described by an oscillation 
around the LO fixed point value. The latter can be estimated from the 
solution of the LO gap equation as described in Sect. |7[l|, which yields 

l.lOMj^f^jrp. Comparing this with the average value for M^(x) for the time 
interval < x° < 60/Minit we find a similar value M = l.lSMpj^jrp. 

One may ask if the agreement can be further improved for not too late 
times by comparing it with an "improved" LO (Hartree) approximation, 
LO"*", that takes into account the local part of the NLO self energy 
contribution (cf. Sect. |5.1| ). In this approximation one neglects the memory 
integrals on the RHS of the NLO evolution equations ( ^.211) and ( ^.22| ). The 



resulting equations are local in time and have the same structure as the LO 
ones, however, with the LO and NLO contribution to the mass term M'^{x) 
included as given by Eq. ( |5.24| ). The large-time limit of the mass term in 
the LO^ approximation is determined by the LO^ fixed point solution in 
complete analogy to the discussion in Sect. |7.1| . The fixed point value can be 
estimated from the corresponding gap equation which yields the asymptotic 
LO"*" value l.llMj^j^jrp. The value slightly improves the LO result. 

Figs. |9HTl] gave an example for the nonequilibrium evolution after a 



"quench" for weak effective coupling where the damping rate is much smaller 

(damp) //- /-;r-r2 



than the frequency of oscillations, i.e. 1 ^ 7q /(VM /27t) = 0.09 for 
the employed coupling \/6N = 0.083 Mj^j^it- In the following we consider a 
"quench" with a larger drop in the effective mass term M^/M^^^^^ = 2.91 
and a stronger effective coupling X/6N = O.lTMj^j^jrp for N = 4. The 
initial particle number distribution is no{p) = l/(exp[-\/p2 _|_ Mq/Tq] — 1) 
with To = 8.47MiNiT- As a consequence of the larger effective coupling the 
correlation functions exhibit stronger damping. The damping rate obtained 
from an exponential fit to the asymptotic behavior of F{x^, 0;p = 0) is found 
to be 'y^^^'^^ = O.llMiNiT- We observe that the oscillation frequency of 
F(x°,0;p = 0) quickly stabilizes around I.IOMinit/Stt = O.ISMinit which is 
of the same order than the damping rate. Correspondingly, one finds from 
Fig. Oan equal-time zero mode F{x^,x^]p = 0) which is effectively damped 



out at NLO after a few oscillations. The deviation from the LO result 
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Figure 12: Time dependence of the equal-time zero mode F{x^,x^;p = 0) 
after a "quench" , similar to Fig. ^ but not in the weak coupling regime (here 
27r^(damp)^^^ ~ C(l)). The inset shows the mass term M (x). 



is substantial. From the inset of Fig. |12| one observes after a short initial 
period an almost constant mass term with a value M^(x) ^ LGMij^^r^. The 
deviation of M^(x), which contains the sum over all equal-time modes, from 
the corresponding LO"*" value of about l.SMj^j^jrp is less severe than for the 
correlation modes. 

We note that the relatively large value for M^(x) infered from Fig. |12 
does not characterize well the oscillation frequencies of the correlation zero 
modes. Following Ref. |T3[ a better characterization for the mode frequencies 
can be obtained from an effective mode energy defined as 

.0 „,o. \ h 



d^ odyoF{x^,y^; py 
F{x^y^-p) 

which coincides with the dispersion relation of a free theory for A — 0. After 
a short oscillation period we find an almost constant value eo(a;'') ~ I.IMjnit- 
One observes that e^lli^ is in good agreement with the oscillation frequency 
for F(x°,0;p = 0) given above. A strongly damped behavior may therefore 
be characterized by 

(damp) 

> 1 • (7.4) 
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Figure 13: Nonequilibrium evolution of the effective four-point coupling at 
NLO for the initial conditions of Fig. |1^ (in units of Minit)- The added 
frame shows the large-time behavior for > 26/Minit on a different scale. 

7.2.3 Effective coupling 

The nonlocal NLO four-point function { f).5\j can be cast into an effective local 
coupling by integration over one variable. For spatially homogeneous initial 
conditions we define the time-dependent effective coupling 

±(1- ly+^y^x^y'-.p)). (7.5) 

One observes that the effective coupling is determined by the antisymmetric 
function Ip{x'^,y^;p) given in Eq. (|5.28|) . 

The nonequilibrium evolution of Xes/QN{x^; p = 0) is shown in Fig. |T^ 
for the same initial conditions as for Fig. |12|. The effective coupling drops 
in response to the "quench" and after a short oscillatory period it exhibits 
a slowly changing behavior. We find a weak dependence of Xes/6N{x^]p) 
on momentum. We note that for the employed Gaussian initial conditions 
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Figure 14: The same large-time behavior of the effective four-point couphng 
at NLO as in Fig. O. Shown is a comparison with the exponential 
fit given in the Figure (all in units of Mjnit). One observes for 
> 250/MiNiT a very good approximation by an exponential with the small 
^ 2.2 X 10"^MiNiT for the approach to equilibrium (cf. Sect. [7.4|) . 



rate 7the 



the time evolution starts with the classical vertex ~ A parametrizing the 
action S, with no radiative or thermal corrections included. At large times, 
when the theory starts to thermalize, the four-point vertex approaches its 
renormalized equilibrium value as is discussed below0. 



7.3 Non-exponential drifting at intermediate times 

We emphasize that the characteristic time scale for damping r^'^^'^p) 
1 / does not correspond in general to the time scale for thermalization. 
Following a "quench" the initial rapid oscillations are effectively damped 
after t^'^^^p), however, the system is typically still far away from equilibrium. 



As an example, the right frame of Fig. 13 shows the asymptotic late-time 



behavior of the four-point function zero mode X^six = 0)/QN for the 

^'^ In this respect the time evolution is reminiscent of the Wilsonian renormahzation 
group flow of couphngs, where the flow interpolates between the classical or microscopic 
parameters and the macroscopic ones p3| . 
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"quench" described above. Its asymptotic time dependence is very well 
approximated by an exponential behavior with rate 7thcrm — 2-2 x 10~'^Minit 
as is demonstrated in Fig. [T^for a wide range of times. This rate characterizes 
the late-time approach to thermal equilibrium (cf. Sect. \7.4\ ). It is much 
smaller than the damping rate and for the employed initial condition the 
corresponding time scales differ by almost two orders of magnitude. 

We also note that the evolution of X^s in Fig. |l^, after time- averaging 
over the oscillation time ~ (^o)"^, cannot be approximated by a simple 
exponential before times of 0(5O/M[nit)- The time dependence is 
parametrically much slower than an exponential and for x° > l/'ylj^^^^^ 
better approximated by a power law behavior. Indeed, for a large variety of 
initial conditions we find a slow change or drifting of equal-time correlation 
modes over time scales Tdnft which are typically much larger than t^^^^p) , 
In the following we address the question in more detail of what happens 
at intermediate times t^^^^^p) ^ 2;° < ^(therm)^ between the two time scales 
^(damp) g^j^^ ^(therm) -y^^j^jcii are both wcU characterized by an exponential rate. 



7.3.1 "Tsunami" 

More than for the above "quench", the characteristic drifting time scale 
becomes apparent for a "tsunami" initial condition (cf. Sect. p. In 
the following we consider an initial particle number no(p) distributed by a 
Gaussian 

no{p) = Aexp (^-^i\p\ - Pts)^^ (7.6) 

peaked around \p\ = pts = SMjnit with a width determined by 
a = O.SMiNiT and amplitude A = 10. The initial symmetric two- 
point function is F(x°,?/°;p)|^o=j^o=o = (^o(p) + 1/2)/ ^/f^TM^O) with 
dxoF{x^,y^;p)\^o=yO=o = and F(x°, 2/°;p) 9^o9yoF(x°, ?/°;p)Uo=j,o=o = 
[no(p) + 1/2]^. The distribution is symmetric under p —p and reminiscent 
of two colliding wave packets in one spatial dimension. The units are 
given in terms of the renormalized initial vacuum mass ( p.24| ) determined 



by MiNiT = M{x)\xO=o for no{p) = 0. Note that this is different from the 
initial effective mass M{x)\x0^q for the nonvanishing particle distribution, 
77-0 (p) 7^ 0, as used above. 



Similar to what is observed for the "quench" in Fig. 0, the unequal-time 
correlation modes F{x^]0;p) and p{x^]0;p) for the "tsunami" in the NLO 
approximation quickly approach an exponential damping behavior. Here the 
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Figure 15: Comparison of the LO and NLO time dependence of the equal- 
time correlation modes F{x^,x^\p) for the "tsunami" initial condition. The 
importance of scattering included in the NLO approximation is apparent: 
the correlation modes drift away from the LO result and the 

"tsunami" decays, approaching thermal equilibrium at large times. 



damping rate for the zero modes is rj^'^™"^^ = 0.022M[nit. If we compare this 
with their oscillation frequency or eo, as discussed below Eq. (|7.3|), one finds 



27r7r"''Veo ^ 0.01 . (7.7) 

The considered "tsunami" is clearly in the weakly damped regime, which is 
due to a small effective coupling A/6A^ = 0.1Mj^j,^it with N = 10. 

In Fig. |15| we present the time evolution of the equal-time correlation 
modes for different momenta. Shown are the LO and NLO 

results for zero momentum, close to the "tsunami" momentum pts and about 
twice pts- One observes that the LO equal-time correlations are constant in 
time. This behavior can be understood from the fact that for the employed 
"tsunami" initial condition the evolution starts at a time translation invariant 



fixed point solution of the LO equations described in Sect. fTl . 

Scattering effects included in the NLO approximation drive the system 
away from the LO fixed point. The equal-time correlation modes can be 
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Figure 16: Effective mode energy ep(x°), defined in Eq. (|7.3|) , for the 
"tsunami" initial condition of Fig. |15[ Comparison with the latter figure 
shows that despite substantial changes in the effective particle number the 
dispersion relation remains almost unchanged. 



interpreted as the ratio of effective particle number^ |]T3[ 

n{x^- p) + ^= [f{x\ y'- p) d,odyoF{x', y'; p)) ' Uo=,o (7.8) 

and effective mode energy ep{x^) as defined in Eq. (|7.3| ), i.e. F(a;°,a;°;p) = 
[n{x^]p) + l/2]/ep{x^). From Fig. ^ one observes that the zero momentum 
mode (dashed line) gets more and more populated as time proceeds, while the 
high momentum mode (dotted line) exhibits relatively small changes. The 
highly populated modes near |pts| become less populated and the "tsunami" 
decays with time (solid line). 

For comparison we show in Fig. 16 the time dependence of ep(x°) for the 
corresponding momenta. As pointed out in Sect. [7.2.2| this quantity gives 
a good description for the oscillation frequencies of the two-point function 
modes. One observes from the figure that after a short oscillatory period and 
a small increase of ep{x^) for low momenta, the effective mode energy is almost 



^^For time-translation invariant two-point functions this definition corresponds to (7.2) 
with dxoF{x" - y°;p)\xO=yO = 0. 
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Figure 17: Double logarithmic plot for the equal-time mode p) close 

to the "tsunami" momentum pts- The dotted straight line corresponds to a 



power law behavior ~ (x^ Mjnit)^'^^- The fit describes F(xO,a;°;p) very well 
for x° < 30/MiNiT if time-averaged over the oscillation period. At later 
times the power law behavior changes to an exponential approach to thermal 
equilibrium. 



constant in time. This is in sharp contrast to the substantial changes in the 
effective particle number as can be infered from Fig. IS. A similar observation 
has been pointed out in Ref. 



13 1 using the three-loop approximation. 



7.3.2 Power law drifting 

We point out that the decay of the "tsunami" modes is not governed by an 
exponential time dependence, here on a time scale r^*^'''^'^ <^ ^(therm)_ rjj^^ 
deviation of the equal-time correlation modes from the LO fixed point is 
rather well approximated by a power law behavior and thus parametrically 
slower than the exponential damping or late-time behavior. In particular, 
we find that the symmetric two-point function behaves approximately as 



(7.9) 



This is demonstrated in the double logarithmic plot of Fig. ^ for 
F{x^,x^]p = 0.94pts)- The NLO curve (solid line) oscillates around, and 
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quickly converges to, the (dotted) straight hne fit characterizing a power 
law. Time-averaged over the oscillation period the evolution of the equal- 
time mode is well approximated by ~ Mjnit)^'''^ for < 30/Minit- 
The details of the intermediate-time power law behavior of the equal-time 
correlation modes depend on the mode momenta. 

The observed behavior is reminiscent of a hydrodynamic regime 
connecting the exponential early-time and late-time regimes. One infers 



from Fig. that at later times a strong deviation from the straight line 
occurs. For times around x° ~ 100/Minit the evolution is governed by 
neither a simple power nor exponential dependence. At sufficiently large 
times it approaches an exponential evolution towards thermal equilibrium 
(cf. Sect. ^). 

7.3.3 Strong interactions at intermediate times 

For higher n-point functions one can observe approximate power law behavior 
for certain intermediate time periods as well, however, the time evolution 
before the exponential late-time behavior can be rather complex. In 
Fig. ^ we show the effective four-point function Aeff(x°;p = 0) for the 



above "tsunami" initial condition. We find that before times of about 
x° < 300/MiNiT the drifting behavior is clearly non-exponential. For later 
times the evolution of AefT(x°;j9) approaches an exponential time dependence 
with a characteristic zero mode rate ^j^*^''''™^ ~ 1.7 x 10~^Minit- The 
thermalization time exceeds the characteristic damping time by more than an 
order of magnitude, similar to what has been observed above for the "quench" 
initial condition. This clear separation of scales happens despite the fact that 



the considered "tsunami" is weakly damped according to Eq. ( |7.7|) 



The effective coupling can, however, change substantially during the 



nonequilibrium evolution. From Fig. |I8] we observe a large increase of 
Aefr(a;°; p = 0) in response to the "tsunami" initial condition. The system 
becomes effectively much stronger interacting for some time before it relaxes 
to equilibrium. For the considered initial condition the phenomenon of a 
strong interaction at intermediate times is restricted to the low momentum 
modes with p ~ 0. For higher momenta we find that the effective coupling 
changes much less between the initial and the late-time value. One observes 
that the quick population of low momentum modes encountered in Fig. 
goes along with a large effective coupling. Any estimate of the involved time 
scale needs to take into account the strong renormalization of the four-point 
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Figure 18: The effective four-point function Xes{x'^;p)/6N at zero momentum 
for tlie "tsunami" employed in Fig. |I5[ The effective couphng for the 
zero momentum mode rises substantially at intermediate times before it 
approaches thermal equilibrium. For x° > 300/Minit the curve is very 
well approximated by an exponential behavior with rate ^^^^"^^^^ ~ 1.7 x 

IQ-^MiNIT- 



function during the nonequilibrium evolution. 

We stress that a simple exponential relaxation can never exhibit such a 
nontrivial behavior. The non-exponential drifting at intermediate times is an 
important phenomenon that separates the early-time and the late-time scales 
associated to exponential damping on the one hand and thermalization on the 
other hand. The overall time needed to reach thermal equilibrium crucially 
depends on this phenomenon and cannot be separated from it. 



7.4 Late-time exponential thermalization 

Independent of the detailed "quench" or "tsunami" initial conditions we 
find for sufficiently large times an exponential behavior to very good 
approximation. This is exemplified for the effective four-point function 
displayed in Figs. |13| and |l^ after a "quench" . The latter figure verifies that 
after a certain time the correlation function follows closely the exponential 
fit for a large time range Ax° = 1000/Minit- A similar observation holds for 
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Figure 19: Effective particle number distribution for a "tsunami" in presence 
of a tliermal background. The solid line shows the initial distribution which 
for low [p <^ pts) and for high {p ^ pts) momenta follows a Bose-Einstein 
distribution, i.e. ln[H-l/n(0; ep)] ~ ep(0)/To with Tq = 4Minit. At late times 
the nonthermal distribution equilibrates and approaches a straight line with 
inverse slope Teq = 4.7Minit > Tq. 



the "tsunami" employed in Fig. |T8|. In particular, for both types of initial 
conditions we find that the corresponding late-time scale exceeds the damping 
time scale by far, r^^^^^^^ ^ ^(damp)_ other words, we observe a very 
efficient damping even for relatively weak coupling or, correspondingly, a 
slow thermalization. 

To observe thermalization on a shorter time scale we discuss in the 
following a "tsunami" initial condition for a stronger coupling X/QN = 
O.SMi^NiT (with Pts = 2.5MiNiT, cr = 0.25Minit and ^ = 4 in Eq. (|T3)). 
In addition, we consider this peaked particle distribution in presence of a 
thermal background with initial "temperature" Tq = 4Minit-|3 The initial 
particle number distribution can be infered from the solid line of Fig. [19 . 
Shown is the combination ln[l + l/n{x^;p)] as a function of ep(x°), where 
the effective particle number and mode energy are defined in Eqs. (|7.8|) and 
( |7.tj| ) respectively. If n{ep) strictly follows a Bose-Einstein distribution with 

i^Thc initial mass term given by Eq. ( |5.24| ) is M(0) = 2.24Minit- 
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temperature T then ln[l + l/n{x'^;p)] = tp/T. Correspondingly, from the 
sohd hne (x° = 0) in Fig. |19| one observes the initial "thermal background" 
as a straight line distorted by the nonthermal "tsunami" peak. In the 
interacting theory for > the employed definitions for effective particle 
number and mode energy should allow for a characterization of the late-time 
thermal equilibrium in a similar way if a "quasi-particle" picture applies (cf. 
also the discussion in Ref. 0). 

Similar to the discussion of the initially weakly coupled "tsunami" in 
Sect. |7.3| , one observes from Fig. |l^ that the initial high occupation number 
in a small momentum range decays. More and more low momentum modes 
get excited and the particle distribution approaches a thermal shape. After 
rapid changes in n(x°; e^) at early times the curves representing snapshots 
at equidistant time steps Aa;° = 10/Minit converge to a straight line to 
high accuracy with inverse slope Teq = 4.7Minit > Tq. The asymptotic 
curve therefore corresponds precisely to what one expects of a thermally 
equilibrated "quasi-particle" distribution. We find that the asymptotic late- 
time result for the distribution is determined by the initial energy density and 
is insensitive to the details of the initial conditions ("thermal fixed point"). 
The latter observation has been explicitly demonstrated before in the three- 
loop approximation of the 2PJ effective action in Ref. |]12| . 

We emphasize again that thermalization or independence of the initial 
conditions cannot be observed in a strict sense because of time-reflection 
invariance (thermal equilibrium is invariant under time translation and 
bears no information about initial conditions). However, our results show 
that thermal equilibrium can be approached very closely without deviating 
from it over times of interest. Of course, the initial conditions for which 
thermalization can be observed in this sense have to comply with standard 
clustering requirements or not too large initial spread in energyQ. It is 
important to note that starting from "first principles" in a time-reflection 
invariant, energy conserving description thermalization can only be achieved 
in the above sense: thermal equilibrium can be approached arbitrarily closely 
with time, indistinguishable for most practical purposes, however without 
reaching it on a fundamental level. This is a matter of principle and not a 
question of the employed approximation. 

^^Note that the energy and its moments are conserved quantities during the 
nonequihbrium time evolution. See also the detailed discussion in Ref. for classical 
field theories. 
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Figure 20: The time dependent inverse slope T{x^]p) of the distribution 
function as shown in Fig. The function T{x^;p) is plotted for low 
momentum po — and high momentum p2 — 2pts, where the slope of 
ln(l + Cp)) is insensitive to small changes in at fixed x° and a 

characterization in terms of a "mode temperature" may be useful. At late 
times T(a;°;p) becomes independent of momenta and the system can be 
described by a global value Teq. 

7.4.1 Inverse slope parameter 

In the following we consider the inverse slope of the effective particle number 
distribution as plotted in Fig. |19|. For low momenta p -C Pts or high momenta 
p ^ Pts, where the slope of ln(l + l/n(a:°; e^)) is insensitive to small changes in 
ep at fixed a characterization in terms of a "mode temperature" T{x^;p) 
may be useful — bearing in mind that away from thermal equilibrium this 
quantity has a precise meaning only in terms of the correlation functions 
employed in the definitions ( [7.3| ) and ( |7.8| ). 

Fig. |2y shows the inverse slope T{x^;p) as a function of time for two 
momenta, po — and p2 — 2pts at NLO (solid line) and LO (dotted line). 
As pointed out above the particle number is conserved at LO and the initial 
nonthermal distribution is "frozen in" . Once scattering is taken into account 
at NLO the distribution evolves. One observes that for the low momentum 
modes T(x°;p) shoots up in response to the "tsunami" initial condition, 
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before dropping to the equilibrium value Teq = 4.7Minit- While the low 
momentum modes get strongly "overheated" at intermediate times with 
respect to the thermal equilibrium value, one finds that the high momentum 
modes exhibit a small initial drop in T[x^]p) which then smoothly evolves. 
The population of the low lying modes involves a rapid change in the particle 
number distribution in which the inverse slope parameter overshoots the late- 
time value. Asymptotically the different "mode temperatures" converge. 

It is interesting to observe that though the low momentum modes deviate 
at intermediate times much stronger from Teq than the high momentum ones, 
there is no substantial difference for the equilibration time as measured e.g. 
by the time at which they have approached Teq within 10 %. A similar 
observation can be made for the "quench" initial conditions discussed above 
where we do not see that some modes thermalize first. 



8 Conclusion and outlook 

The generating functional for 2PI Green's functions provides a powerful 
technique for out-of-equilibrium physics. In combination with a 
expansion it yields a controlled nonperturbative description from "first 
principles" which is time-reversal invariant and energy conserving. The 
approach is not restricted to weak couplings or situations close to thermal 
equilibrium or to an effective description based on a separation of scales. 
With the present approach we have a quantitative method at hand which 
can deal with far-from-equilibrium processes, where the latter assumptions 
are typically not justified. An example is the possible strong renormalization 
of the effective four-point function at intermediate times for weak initial 
coupling (cf. Sect. [7.3.3| ), which stresses the fact that a clear separation of 
scales can be difficult to achieve away from equilibrium. 

In this work we have studied two classes of initial conditions representing 
a "quench" and a "tsunami", both in the weak (27r7Q*^^™^^/eo -C 1) and 
the stronger coupling regime 

(27r7j'^'''"PVeo - 0{1)). Independent of the 
details of the initial conditions we observe three qualitatively different time 
regimes, as summarized in Fig. In particular, we find t^'^^^p) <^ ^(therm) 
which emphasizes that the damping time cannot be identified with the 
thermalization time. It is possible to relate the damping rate to the "width" 
of the Wigner transformed spectral function p in the same way as discussed in 



Ref. [O], however, the thermalization rate cannot. Processes with nontrivial 
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momentum exchange are necessary for thermalization and are taken into 
account at NLO which includes off-shell effects^. On-shell particle number 
changing diagrams as the one in Fig. ^ which appears at next-to-next-to- 
leading order, can be important since they contain dynamics on which the 
bulk viscosity depends p|, |29|. 

We emphasize that scattering effects are crucial, even for rather early 
times. None of the features summarized in Fig. |l] are captured by the 
mean field type LO approximation which neglects direct scattering; no 
exponential damping (in particular, no suppression of correlations F{x^, 0;p) 
and p{x^,0;p) with the initial time), no drifting and no thermalization can 
be observed in that case. The LO approximation typically breaks down after 
the first few oscillations and is approximately valid only for x° <^ ^(damp) 
However, for suitably time- or momentum-averaged quantities and not too 
strong couplings we find that the LO fixed point solution provides a rather 
accurate estimate up to times ~ ^(damp) rpj^^ latter observation agrees 
with what has been found for the corresponding classical field theory in 



Ref. 32 



We note that we are able to observe the same qualitative properties 
using the three-loop approximation of the 2PI effective action for = 1 
as employed in Refs. [|^, . Despite quantitative differences this points 
out that the quahtative properties discussed here do not depend on the 
detailed implementation of scattering and memory effects. In particular, 
this emphasizes the importance of the setting sun type diagram of Fig. |^. A 
detailed quantitative comparison is deferred to a separate publication which 
also discusses universality (independence of initial conditions) and parametric 
dependences of time scales pO| . 

There are many directions to pursue which are beyond the scope of the 
present work. With a more efficient numerics — the computations were done 
on a PC with unparallelized code — one can directly address similar problems 
in higher dimensions including the possibility of spontaneous symmetry 
breaking. This would allow for a quantitative quantum field theoretical 
study of the formation of disoriented chiral condensates in the context of 
heavy-ion collisions or of (p)reheating at the end of inflation in the early 
universe along the lines presented here. In the present form the time- 

^^In 1 + 1 dimensions on-shell two-to-two scattering is constrained by the energy 
conservation relation such that it does not change the particle numbers for the involved 
momentum modes. 
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nonlocal scattering terms are difficult to estimate analytically. Though their 
numerical implementation is rather simple the calculations are time and 
memory consuming. However, one may expect that in higher dimensions 
both damping and thermalization become more efficient which reduces the 
need to follow the memory kernels to very large times. It is important to 
have a controlled, quantitative method which can serve as the basis for the 
development of approximation schemes which may allow to find reliable 
analytical solutions in limiting cases. The present techniques present a 
wealth of possibilities to study the applicability of standard phenomenological 
approaches. The range of validity of the classical Boltzmann equation 
or more sophisticated kinetic descriptions can be tested as well as other 
typical approximation schemes based on local equilibrium or linear response 
assumptions. 

We emphasize that the current approach is not restricted to scalar 
fields though the presence of a small nonperturbative expansion parameter 
depends, of course, on the respective theory. It would be very desirable to 
extend the approach to include fermions and gauge fields. 

Acknowledgements 

I thank G. Aarts and J. Cox for collaboration on related work, and 
A. Jakovac, T. Prokopec and C. Wetterich for interesting discussions. 
I would hke to thank W. Wetzel for support with computers. 

References 

[1] J. Schwinger, J. Math. Phys. 2 (1961) 407; L. V. Keldysh, Zh. Eksp. 
Teor. Fiz. 47 (1964) 1515 [Sov. Phys. JETP 20 (1965) 1018]. 

[2] For a digrammatic analysis in the context of a scalar theory see e.g. 
S. Jeon, Phys. Rev. D52 (1995) 3591; S. Jeon, L. Yaffe, Phys. Rev. 
D53 (1996) 5799. 

[3] For a review see D. Bodeker, Nucl. Phys. Proc. Suppl. 94 (2001) 61, and 
references therein. 

[4] Among the extensive literature on the leading-order approximation 
see e.g. F. Cooper, S. Habib, Y. Kluger, E. Mottola, J.P. Paz, P.R. 
Anderson, Phys. Rev. D50 (1994) 2848; D. Boyanovsky, H.J. de Vega, 



46 



R. Holman, D.S. Lee, A. Singh, Phys. Rev. D51 (1995) 4419; S. Habib, 
Y. Kluger, E. Mottola, J.R Paz, Phys. Rev. Lett. 76 (1996) 4660; 
D. Boyanovsky, H.J. de Vega, R. Holman, J.F.J. Salgado, Phys. Rev. 
D54 (1996) 7570; F. Cooper, S. Habib, Y. Kluger, E. Mottola, Phys. 
Rev. D55 (1997) 6471; J. Baacke, K. Heitmann, C. Patzold, Phys. 
Rev. D55 (1997) 2320; Phys. Rev. D57 (1998) 6406; D. Boyanovsky, 
C. Destri, H.J. de Vega, R. Holman, J. Salgado, Phys. Rev. D57 (1998) 
7388; D. Boyanovsky, H. J. de Vega, R. Holman, J. Salgado, Phys. Rev. 
D59 (1999) 125009; C. Destri, E. Manfredini, Phys. Rev. D62 (2000) 
025007; and references therein. 

[5] G. Aarts, J. Smit, Phys. Rev. D61 (2000) 025002; M. Salle, J. Smit, 
J. C. Vink, |hep-ph/001236^ , |hep-ph/0012346| . L. M. Bettencourt, 
K. Pao, J. G. Sanderson, |hep-ph/01042TO| . 



[6] 
[7] 



[9 
[10 

[11 

[12 
[13 
[14 

[15 
[16 



L.M.A. Bettencourt, C. Wetterich, Phys. Lett. B430 (1998) 140; 

B. Mihaila, J. F. Dawson, F. Cooper, Phys. Rev. D56 (1997) 5400; 
L. M. Bettencourt, C. Wetterich, |hep-ph/9805360| ; B. Mihaila, T. Athan, 

F. Cooper, J. Dawson, S. Habib, Phys. Rev. D62 (2000) 125015; 
A. V. Ryzhov, L. G. Yaffe, Phys. Rev. D62 (2000) 125003. 

J. M. Cornwall, R. Jackiw, E. Tomboulis, Phys. Rev. DIO (1974) 2428. 

K. Chou, Z. Su, B. Hao, L. Yu, Phys. Rept. 118 (1985) 1. 

E. Calzetta, B. L. Hu, Phys. Rev. D37 (1988) 2878. 

See also J.M. Luttinger, J.C. Ward, Phys. Rev. 118 (1960) 1417; 

G. Baym, Phys. Rev. 127 (1962) 1391; B. Vanderheyden, G. Baym, 
|hep-ph/0002"29ll . 



J. Berges, J. Cox, |hep-ph/0006T60 



G. Aarts, J. Berges, |hep-ph/0l030ig| . 

L.P. Kadanoff, G. Baym, "Quantum Statistical Mechanics", Benjamin, 
New York (1962). 

P. Danielewicz, Ann. Phys. 152 (1984) 239; 197 (1990) 154; 
S. Mrowczynski, P. Danielewicz, Nucl. Phys. B342 (1990) 345; 



47 



[17] S. Mrowczynski, U. Heinz, Annals Phys. 229 (1994) 1. 
[18] J. Blaizot, E. lancu, |hep-ph/OlOlTUB| . 

[19] Y. B. Ivanov, J. Knoll, D. N. Voskresensky, Nucl. Phys. A672 (2000) 
313; S. Leupold, Nucl. Phys. A672 (2000) 475. 

[20] M. Joyce, K. Kainulainen, T. Prokopec, Phys. Lett. B474 (2000) 402. 

[21] J. Herges, K. Rajagopal, Nucl. Phys. B538 (1999) 215; M. A. Halasz, 
A. D. Jackson, R. E. Shrock, M. A. Stephanov, J. J. Verbaarschot, Phys. 
Rev. D58 (1998) 096007. 

[22] H. Mihaila, F. Cooper, J. F. Dawson, Phys. Rev. D63 (2001) 096003. 

[23] K. Rajagopal, F. Wilczek, Nucl. Phys. B404 (1993) 577. 

[24] R.D. Pisarski, |hep-ph/9710370| ; D. Hoyanovsky, H.J. de Vega, 
R. Holman, S.P. Kumar, R.D. Pisarski, Phys. Rev. D57 (1998) 3653. 

[25] L. McLerran, R. Venugopalan, Phys. Rev. D49 (1994) 2233; (1994) 
3352; for a recent review see e.g. L. McLerran, |hep-ph/0104285| . 

[26] R.D. Pisarski, private communication. 

[27] P. Arnold, L. G. Yaffe, Phys. Rev. D57 (1998) 1178. 

[28] F. Cooper, S. Habib, Y. Kluger, E. Mottola, Phys. Rev. D55 (1997) 
6471. 

[29] E.A. Calzetta, B.L. Hu, S.A. Ramsey, Phys. Rev. D61 (2000) 125013. 
[30] G. Aarts, J. Herges, in preparation. 

[31] C. Wetterich, Phys. Rev. E56 (1997) 2687 ; Phys. Rev. Lett. 78 (1997) 
3598. 

[32] G. Aarts, G.F. Honini, C. Wetterich, Phys. Rev. D63 (2001) 025012. 

[33] For a recent review on this subject see J. Herges, N. Tetradis, 
C. Wetterich, appears in Phys. Rept. ( |hep-ph/0005122| ) . 

[34] G. Aarts, G. F. Honini, C. Wetterich, Nucl. Phys. B587 (2000) 403 . 



48 



